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We prove stability of the spectral gap for gapped, frustration-free Hamiltonians under general, 
quasi-local perturbations. We present a necessary and sufficient condition for stability, which we 
call Local Topological Quantum Order and show that this condition implies an area law for the 
entanglement entropy of the groundstate subspace. This result extends previous work by Bravyi 
et al. on the stability of topological quantum order for Hamiltonians composed of commuting 
projections with a common zero-energy subspace. We conclude with a list of open problems relevant 
to spectral gaps and topological quantum order. 



I. INTRODUCTION 



Recent interest in topological quantum computation has focused the attention of the condensed matter and math- 
i— i ematical physics community on Hamiltonians whose low-energy sectors exhibit some form of topological order. In 
a seminal paper by Kitaev [25] . a Hamiltonian known as the toric code, was constructed out of commuting spin- 
interaction terms, such that the groundstate subspace exhibits four-fold degeneracy on a torus, which may only be 
detected under macroscopic operations. In particular, the different groundstates are indistinguishable on microscopic 
and even macroscopic scales, which implies that local errors cannot create logical errors in any encoding that utilizes 
such topologically ordered groundstates as the encoded state of two qubits. 

Nevertheless, if under a local perturbation the groundstate subspace were to split appreciably or become mixed 
1—1 with the non-topological higher-energy sectors, then any guarantee of protection from errors would no longer be valid. 
At this point, another type of stability against errors is required, this time at the level of the spacing between the 
low-energy and high-energy sectors and in the form of a negligible splitting in the groundstate subspace. Progress 
in this direction was made by Klich |26j using the method of cluster expansions, and independently by Bravyi et 
aZ. [51 IE] ■ In particular, the latter result showed that all Hamiltonians with commuting spin-interaction terms have 
negligible splitting in the groundstate subspace and a spectral gap that is stable under local perturbations, as long as 
the groundstates satisfy some type of local indistinguishability and frustration- freeness condition. Since the toric code 
Hamiltonian satisfies both conditions for gap stability, it follows that encoded qubits in the toric code groundstate 
subspace are robust against local errors, even if the Hamiltonian interactions are not precisely engineered. Moreover, 
at low enough temperatures, the spectral gap implies that the creation of excitations is unlikely to happen in relatively 
small time scales. 

Even so, the toric code has a seemingly fatal flaw: at non-zero temperature, local excitations, once created, can 
travel around the torus to create logical errors at no extra energy cost. Motivated by this issue, soon after the results of 
Klich and Bravyi et al., a series of papers appeared that focused on the beneficial effects of Anderson- type localization 
of the undesired excitations in the presence of impurities, or external magnetic fields [Ml ESI E9 • 

Still, until recently, no Hamiltonian model was known that combined the desirable properties of the toric code, with 
a rigorous and sufficiently large lower bound on the energy barrier restricting the mobility of unwanted excitations 
at non-zero temperature. In fact, it was shown by Yoshida [46], that such Hamiltonians would need to either forgo 
commutativity or break certain "natural" conditions, such as translational invariance or scale invariance, in order 
to satisfy stability at non-zero temperature. Indeed, the only known family of Hamiltonians with all the desired 
properties, recently presented by Haah in |14j . has logical operators with fractal geometry in three dimensions and no 
scale invariance, thus sidestepping Yoshida's no-go theorem. Yet, despite the rigorous bound [5] on the energy barrier 
of Haah's Hamiltonian, there remains a question of whether the barrier is large enough to allow for operations on the 
logical qubits, like read-out and error-correction, which may require times comparable to the time it takes for logical 
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errors to appear. Moreover, it was recently shown by Hastings [IB], that all two-dimensional Hamiltonians which are 
a sum of commuting terms, have no topological order at non-zero temperature. 

Motivated by this line of research and the larger question of the classification of quantum phases [31 [SJ [3j|l SI] > 
we present here a generalization of the result by Bravyi et al, which removes the commutativity of the Hamiltonian 
terms as an assumption for stability. Some of the new candidate Hamiltonians now include, parent Hamiltonians 
of Matrix Product States (MPS) and Projected Entangled Pair States (PEPS) [T3J [36l HQ1 EI] , as well as all other 
frustration-free Hamiltonians - that is, Hamiltonians whose groundstates minimize the energy of each local interaction 
term. Moreover, we generalize the conditions needed for the stability of the spectral gap, in hopes that in the future, 
one may be able to prove an equivalent result for general, gapped Hamiltonians, whose groundstates satisfy some type 
of topological order. 

The stability of quantum phases was already studied by Borgs et al. [1| and Datta et al. [TUl [H] , where "classical" 
systems were shown to be robust against small quantum perturbations up to some low-temperature, using the methods 
of contour and cluster expansions. Nevertheless, here, we draw heavily from the methods developed in [6,8, following 
the more succinct format of [6], which uses Hastings' powerful quasi- adiabatic continuation [21]. We make extensive use 
of Lieb- Robinson bounds [57] , both for the evolution of operators according to the Hamiltonians under consideration 
[TSJ, [T13 , 29 33; (37J [35], as well as Lieb- Robinson bounds on the quasi-adiabatic evolution, or spectral flow, of gapped 
eigenspaces [31 H71 [20] , 

The paper is organized as follows: In Sec. [H] we define the class of Hamiltonians whose stability we proceed to 
study. In Sec. |III[ we introduce the conditions sufficient for proving stability, clarifying with examples the extend 



to which the conditions are also necessary. In Sec. |IV[ we define the class of perturbations we allow and proceed to 
transform the perturbed Hamiltonian, through a unitary process and a global energy shift, in a form amenable to 
studying its low energy sectors. Sec. |V] contains the proof of the bound on the energy of the perturbation with respect 
to arbitrary states, relative to the energy of the unperturbed Hamiltonian. The main theorem appears in Sec. [VI] and 
in Sec. VII we close with a discussion of future work and a list of open problems. Finally, Appendix |A"] contains the 
technical details behind the estimate in Lemma [21 



II. THE HAMILTONIAN 



We study the stability of the spectral gap for Hamiltonians Hq defined on A = [OjL]'' C Z d , satisfying: 

1. (Spatially-local) Hq — Y2 u eA Qu, where each interaction Q u is a projection acting non-trivially on 6„(1), where 
b u (r) is the ball of radius r with center at site u £ A. 

2. (Periodic-boundary) The set A satisfies periodic boundary conditions. 

3. (Frustration-free) For Pq the projector onto the groundstate subspace of Hq, we have HqPq = 0. 

4. (Gapped) H Q has a spectral gap 7l > 7 > 0, for all L>2. 

Recall that L is the linear size of the support of Hq, so when we write Hq, we assume a fixed size L. The gap 
condition implies that as Hq is defined (translationally invariant, or otherwise) for larger L, the gap remains 
uniformly bounded from below. The restriction of each interaction on a ball of radius I is not crucial, due to coarse- 
graining, as long as the interactions have constant range of support. Moreover, the assumption that the interactions 
are projections is not crucial either, as long as Q u Pq — X u Pq, where X u is the minimum eigenvalue of Q u . This 
condition is the most general way of defining a frustration-free Hamiltonian. From now on, however, we will assume 
that X u = 0, which follows from substituting Q u with Q u — X u l. Notice that this transformation corresponds to a 
global energy shift by — y] X u 1, which leaves the spectral gap, and all other relevant properties of the Hamiltonian, 
unchanged. 

III. ASSUMPTIONS FOR STABILITY 

We begin this section by defining the following projections onto low-energy eigenstates: 

Definition 1. For B — b u {r), r > 2, we define i"s(e) to be the projection onto the subspace of eigenstates of 
Hb '.= Y2b (i)cb Q"" w ith energy at most e > 0. We set Pb '■= -Pb(O). 

The above definition implies the following relations, which hold for all frustration-free systems: 

Corollary 1. The following claims follow from the definition o/Pb(e); 
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• H B >e(l-P B (e)). 

• P B (e)P B (e / ) = P B (e'), for < e' < e and P A (e) = P , /or < e < 7 . 

• P B (e) P c (0) = P c (0) } forBcC and e > 0. 

• P B (e) = Pb(0) if Hb has spectral gap greater than e. 

Proof. The first claim follows directly from the definition of P B (e), since 1 — P B (e) includes only states with energy 
at least e in Hb- The second claim follows from the fact that 1 — P B (e) is orthogonal to Pb(c') for e' < e. Moreover, 
Pv( e ) = Pa(c') = Po, for all e' € [0, e], since Hq has spectral gap 7 > e. The third claim follows from the second claim 
and Pb(0)Pc{0) = Pc(0) ; which is true for frustration-free Hamiltonians since H c Pc{0) — implies H B Pc(0) = 0. 
In particular, we get: P B (e)P c {0) = {P B {e)P B (0)) P c (0) = Pc(0). The fourth claim follows from the definition of 
P B (e) and the spectral decomposition of Hb- □ 

We are now ready to introduce the main assumption needed for the stability of frustration-free systems. It is a 
condition on the degree of topological order for local ground states, so we shall call it Local Topological Quantum 
Order. 



A. Topological Order Assumption. 

Assumption 1. [Local-TQO] Let A — b u (r) with r < L* < L and Oa be any bounded operator supported on A. Set 
A(£) = b u (r + £) and define 

c e (0 A ) := Tr(P AW A )/TrPA(e), <O a ) := Tt{P q O a )/ P . 

Then, for each fixed 1 < £ < L — r, we have: 

\\P A{t) A P A{t) -ci{0 A )P m \\ < \\O A \\A {£), 

where Ao(^) is a decaying function of I. 

The cut-off parameter for topological order, L* , may be thought of as the code distance of the encoding one 
may define on the (degenerate) groundstate space of Hq. Since we are interested in the stability of the spectral 
gap of Hq, we expect that L* should scale with L. To see this, note that if the Local-TQO condition fails after 
some fixed distance, then there must be operators O a , with A of increasing diameter and \\O a \\ = 1, such that 
(^J, Oa^q) — (^q, Oa^o) — c 0; f° r some constant c > and {|^o) > |*o)} groundstates of H . If we assume that 
the number of these operators Oa, for a given size of their support A, scales with L, then it follows that perturbing 
^0 by J2a®a, would create a macroscopic energy gap within the groundstate subspace, thus destroying any notion 
of stability and macroscopic order. This is exactly what happens in the case of the two-dimensional Ising Hamiltonian 
Hq = — X) p ^ 9 a p cr q, under a uniform magnetic field JJ2 P a p- -"- n particular, the two- fold degenerate groundstate space 
fails to satisfy Local-TQO at the most basic level, since any local operator a*, p E Z 2 can distinguish the groundstates 
(one has eigenvalue +1 and the other —1.) 

One may wonder at this point why it is necessary to demand that low-energy states are locally indistinguishable at 
all scales, and not just at the global scale of the system. In other words, why do we care if the groundstates of the local 
Hamiltonians are also topologically ordered? The reason is best exemplified with the following counterexample: Let 
H = — erg — J2 P ~q cr p <T q: where p, q e Z 2 . This is a gapped Hamiltonian with a unique groundstate |0 • • • 0), satisfying 
topological order at the global level. It is easy to see that a perturbation of the form JJ2 P eA a p' wu ^ c l° se the g a P 
between the states 1 000 • • • 0) and 1 1 1 1 ••■!), even for J ~ 1/|A|. Since Local-TQO fails for any region that does not 
include the origin, this implies that some form of Local-TQO is necessary to protect the spectral gap from collapsing. 
The above examples validate our intuition that Hamiltonians are unstable because local "order parameters" can act as 
perturbations to either open the gap between groundstates, or close the gap between global groundstates and excited 
states which have low energy locally. In that sense, the above condition implies that states that minimize the energy 
of the Hamiltonian as more and more interactions are added, differ only near the boundary of the region we are 
considering. 

The Local-TQO condition defined above, combines and generalizes assumptions TQO-1,2 of the stability result for 
the commuting case |6l|8j. It is worth noting that the above condition is a property of the local, low-energy subspaces. 
In fact, it is possible in many cases (e.g. Hamiltonians with a unique Matrix Product State as their groundstate, or 
injective PEPS [3B], @D], [H|) to modify a frustration- free Hamiltonian Hq to a different Hamiltonian Hi, satisfying 
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Local-TQO without changing the global ground state subspace, or closing the gap. In such cases, Hq and Hi are 
both gapped, frustration-free and share the same ground state space, so H s = (1 — s) H + sHi is also gapped and 
frustration-free with the same ground state space. In other words, if Hi satisfies Local-TQO and H does not, we can 
apply the stability result to Hi to infer properties about the groundstate of Hq. For a concrete example, one can look 
at two different toric code Hamiltonians, as discussed in Sec. 2.4]. 

Corollary 2. Let Oa be any bounded operator supported on the ball A — b u (r), with r < L* and u G A. Then, 
Local-TQO implies the following bounds for 1 < £ < L — r: 



\\\O a P aw \\ - \\0 A Po\\\ < 2\\O a \\ y/A {£). (1) 
\\P A (i) - Pa(i)Pa\\ < V3y/A (£). (2) 
Proof. From Local-TQO and the triangle inequality, we have: 

\\\o A P m \\ 2 -c e (\o A \ 2 )\ < \\O A \\ 2 A (£), (3) 

I \\O A P \\ 2 - c L _ r {\0 A \ 2 )\ < \\Oa\? A (L - r), (4) 
\cl-A\Oa\ 2 ) - c e (\0 A \ 2 )\ < \\0 A \\ 2 (A (£) + A (L - r)). (5) 

The last inequality follows from a triangle inequality on Q and the bound: 

\\Po\0 A \ 2 Po ~ c e (\O a \ 2 ) Poll = \\P Pa W [\Oa\ 2 ~ d (\0 A \ 2 )] P A (i)Po\\ < \\PA(t)\0 A \ 2 Pa {1 ) - c e (|Ot| 2 ) P A{t) \\. 

Recalling that A (£) is decaying, we have: 



< 2||O a || 2 (A (£) + A (L - r)) < 4||O a || 2 A (£). 

AW I 



\\\o A p A(l) \\ - \\o A p \\\ 2 < \\O a P aw \\ - \\O A P \\ 2 

The second bound is equivalent to showing that ||(1 — P A )P A (e)\\ 2 < 3A (l). If we are given two projections P and Q 
satisfying: \\PQP — P\\ < S, then we have ||P(1 — <3)P| < <5 => | (1 — Q)P|| 2 < S. In our case, we have from Local- 
TQO: \\P AW P A P A (t)-P A (e)\\ < 2A (i)+A (L-r), where we used Corollary[l]to get c L - r (P A ) = Tr(P A P )/TrP = 1 
and combined it with (|5| and a triangle inequality. □ 

Note that in the above Corollary, the projections P A and P A (i) correspond to Pa(c) and Pa(£)(£), respectively. 
All that is needed for the bounds to hold is that the projections satisfy the Local-TQO condition. Of course, for 
Pa = Pa{0) and P A{e) = Pa(£)(0), we have from Corollary [ll that \\P A (e) - Pa{1)Pa\\ = 0- 

The next Corollary makes precise the notion of local inaistinguishability between ground states and provides an 
equivalent condition to Local-TQO. 

Corollary 3. Let A = b u (r), with r < L* and u E A. For any two groundstates \ijji(A(£))) , \ip2(A(£))) of H A ^; that 
is, P A (i) \ijji(A(£))) = \ipi(A{£))) , i = 1,2, define the reduced density matrices Pi(A) = Ti A m\A |'0i(A(£)))(^(A(^))| . 
Then, Local-TQO is equivalent to: 

\\pi(A)- P2 (A)\\i<2A (£). (6) 

Proof. First, we prove that Local-TQO implies the above condition. Set pi = \ipi(A(£)))(ipi(A(£))\, i = 1,2. Then, 
Pi(A) = r lT: A (t)\A Pi and is supported only on A. From Local-TQO, for any operator Oa supported on A: 

| Tr( Pl - P2 )0 A \ = | (MA(fi)\ °a - ci(Oa) \MMQ)) - (MMQ)\ °a - c e (0 A ) \MMM I < 2\\O A \\A (£). 

Now, recalling that ||pi(A) — /0 2 (^4) 1 1 1 = su Po A -\\o A \\=i I Tr(/°i — P2)O a \, we get the desired bound. To prove the other 
direction, we note that it is sufficient to show that for \ 1 $ (A(£))) satisfying P A (e) I^oC^CO)) = l^o(-^W))i we have: 

| (* (A(£))\O A \^o(A(£))) - c e (0 A )\ < 2||O A ||A (£). (7) 

Setting N = TrP AW -l andp*(A) = Tr m \ A \^ k (A(£)))(^ k (A(£))\, < k < N, where {\V k (A(£))) (V k (A(£))\}% =Q 
is an orthogonal decomposition for Pa(i), we have for ||Oa|| = 1 ; 

\(y (A(£))\O A \* (A(£))) - Ci (0 A )\ = |Tx APo (A) - a{0 A )\ = 

N , N 



J ±^\J2TvO A (po(A)-p k (A))\ < -J—J2\\ PQ (A)-p k (A)\\i < ^AoW- 

k=0 k=l 

Hence, up to constant pre-factors, (|6| and Local-TQO are equivalent. □ 
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Moreover, Local-TQO implies an area law [TB] for the entanglement entropy of the ground states of Hq. 

Corollary 4. Let A — b u (r), with r < L* and u <E A. Any groundstate |\&o) °f Ho satisfying Local-TQO, also 
satisfies an area law for the entanglement entropy of p A ■— Tr^c |^o)(^ r o|- 

S(p A ) <(c d lnD)(l + r) d - 1 -£ , (8) 

where c d is a constant depending only on the dimension d of the lattice, D is the maximum dimension of the on-site 
Hilbert spaces (e.g. D = 2 for spin-1/2 particles) and £q = min{£ : Ao(^) < ^/(l + r)}. 

Proof. We begin by writing the Schmidt decomposition of I'I'o) over the partition of the lattice A(£) : A\ A(£), where 
I will be chosen later: 

l*o) = X)v^*l**(^))>®l**(A\^))>, (* k (A\A(£)),* l (A\A(£))) = 6 k , l , = 1, u k > 0. (9) 

k>0 k>Q 

Note that for frustration-free Hamiltonians, the states {|^ , fc(A(^)))} / t>o in the Schmidt decomposition of \^o) are 
groundstates of the Hamiltonian H A ^y To see this, recall that Hq = 5Z u Qu, Qu > and H \^f ) = 0, which 
implies that (* \H m |tf„) = = E fc >o ^ (* k{A{£))\ H A(e) \V k (A{£))) , hence H A{e) \* k (A{£))) = 0, k > 0. Setting 
p fc (A) = Ti AW \ A \* k (A(t)))(* k (A(t))\, we get from §: 

p A := Tr |* )(*o| = <°*( A ) = ^( A ) + E ^ (^( A ) ~ <°°(^))- (10) 

N fc>0 k>0 

Corollary[3]and (JlOj give: ||^ - Po(^)||i < 2A (€). Combined with the bound S(p (A)) < (\A(£)\ - \A\) hiD < 
d d (l + ry^^£ \nD, for I < r and a constant depending only on d, we have: 

S(ft0 < S( Po (A)) + \S( PA ) - S(p (A))\ < c' d hxD (1 + r) d ^ + eg InD A„(i) (1 + r) d + JT(A (*)), (11) 

where we used the Fannes-Audenaert inequality [2j [12] for the entropy: 

\S( PA ) - S(p (A))\ < (\\ PA - Po (^)|| 1 /2) InD^ + H(\\p A - ft)(A)||i/2), 

with = -plnp - (1 - p) ln(l -p), P e [0, 1]. Setting £ = min{^ : A (f) < + r)} and c d = max{<4 c^'} + 1, 
we get the desired bound, noting that H(A (£ )) < In 2. □ 

We note that for stabilizer Hamiltonians, like the toric code, the Local-TQO decay Ao(£) vanishes after a constant 
length £q; in other words, Aq(£) = 0, I > £q. This implies that for such Hamiltonians, the groundstates satisfy an 
area law exactly. On the other hand, if the decay Aq(£) is exponential in £, then there is a logarithmic divergence 
(£q ~ In (1 + r)) in the area law, which is natural for gapless systems (the above proof made no use of a spectral gap.) 



B. Local Gap Assumption 

The second assumption imposes a condition on the spectral gaps of the local Hamiltonians Hf,t r \. 

Definition 2. [Local-Gap] We say that H Q is locally gapped with gap "f(r) > 0, iff for each u £ A and r > 0, 
Pb u (r){l{ r )) = Pb u (r)- When 7(7*) decays at most polynomially in r, we say that Hq satisfies the Local-Gap condition. 

Below, we present an example of a gapped, translationally-invariant, frustration-free Hamiltonian with Local-TQO, 
that is unstable under vanishingly small perturbations: 

Example 1. Let H 2N = £™ H k . k+1 , where H k . k+1 = l/(3JV)|0)<0| fc <8 |l)<l| fc+ i + |l)<l| fc <8> |l><l| fc+ i + |0)(0| fc ® 
|0)(0|fc + i, fork even and H Kk+1 = 1/(3JV)|1) (l\ k ® |0) (0| fc+1 + 11) (l\ k ® |1) (l| fc+1 + 10) (0| fc <E> |0) (0| fe+1 , for k odd. There 
is a unique groundstate 1 010L . . . 01) , with energy and spectral gap 2/3 to the first excited state 1 1010 . . . 10). Note 
that since all eigenstates of H 2 n can be written as product states over |0) , any state with consecutive 's, or 1 's has 
energy at least 1. It should be obvious that this Hamiltonian satisfies frustration-freeness and Local-TQO (trivially, 
since the groundstate is unique and is a product state), but is unstable under the perturbation V = — Y^j—i Vfc » 
where V k = l/(3A r )|0)(0| fe <g) |l)(l| fc+ i, for k even and V k = l/(3JV)|l)(l| fc <g> |0)(0| fe+ i, for k odd. Interestingly, the 
Hamiltonian H 2 n + V is frustration-free, has degenerate ground state space and satisfies the Local-Gap condition, but 
not the Local-TQO assumption, since there exist local operators, like |0)(0| <g) |1)(1|, that distinguish the two ground 
states {|010101 . . . 01) , 1 101010 . . . 10)}. In other words, the two assumptions for stability seem to be independent and 
dual to each other. Moreover, this example suggests that, in some further refined form, both conditions are necessary 
for stability. 
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Such violations to the Local-Gap condition are not obvious for Hamiltonians that are sums of non- vanishing terms. 
In fact, to the authors' knowledge, it remains an open problem whether it is possible to construct examples of gapped, 
frustration-free Hamiltonians that are sums of local projections and do not satisfy the Local-Gap condition. Moreover, 
we expect that the Local- Gap condition can be merged with a stronger version of the Local- TQO condition, where 
the projections Pa(£) ar e substituted with P^m(e), for some e < 7. Furthermore, we expect that using the powerful 
detectability lemma pQ, one may be able to remove the Local- Gap assumption, when the Hamiltonian interactions can 
be rescaled to projections in a system-size independent way. 

IV. TRANSFORMING THE PERTURBED HAMILTONIAN 

We begin this section by defining the general class of perturbations we will be studying. The only restriction is a 
certain notion of quasi-locality that we make precise in the following definition: 

Definition 3. We define a perturbation V to have strength J and decay fo(r), iff we can write V — X)«eA Zir>o Vv.(j), 
such that V u {r) has support on b u (r) and satisfies ||V^(r)|| < J/o(r), for some rapidly decaying function fo(r). We 
write that V is a (J, fo) -perturbation. 

We have chosen to leave the decay rate of fo ill-defined, beyond saying that it is rapidly-decaying, in order to avoid 
constraining the results that follow unnecessarily. Of course, throughout this paper, we assume that fo decays fast 
enough to allow for the usual Lieb- Robinson bounds [T§1 13"3] on evolutions based on Ho + V. This implies that fixed 
range and exponentially-decaying perturbations are well within the purview of this paper. 

In order to relate the spectrum of the initial Hamiltonian H to the spectrum of the final Hamiltonian H + V, we 
define the following one-parameter family of gapped Hamiltonians: 

Definition 4. Let H s := Hq + sV , where Hq satisfies the assumptions of the previous sections and V is a (J,fo)- 
perturbation. We assume that H s has spectral gap 7 S > 7' > and groundstate projection Po(s), for s G [0, s'] and 
s' G [0, 1] the maximum value of s, such that the gap condition is satisfied. 

Note that s' may depend on 7', J and fo, as well as the linear size L, of the system. Moreover, 7' is fixed a priori 
to some value less than 7. Concretely, we may take 7' =7/2, but this is not crucial. The estimates below will depend 
implicitly on our choice of 7'. 

The main goal of this section is to transform a ( J, / )-perturbation V, into a ( J, w)-perturbation W with similar 
decay, such that: W u (r)Pj, u r r ) — and for some unitary U we have, W (Hq + V)U = Hq + W + A + El, where 
||A|| — > as L — > 00 and E — Tr(UP W (H + V))/TrP . Recall that P&„( r ) is the ground state projector of the 
local, frustration- free Hamiltonian H bu ( r y 

A. Making the global ground states commute quasi-locally. 

In this section, we re-derive a version of Lemma 7 in [6], which will be sufficient for our purposes. The content of 
this Lemma can be summarized as the decomposition of any gapped Hamiltonian into a sum of quasi-local interactions 
that commute with the groundstate subspace. Note that any local, frustration-free Hamiltonian satisfies this condition 
(with strict locality), but once we add a perturbation, the new terms may no longer commute with the groundstate 
subspace. 

Lemma 1. Let H s be the family of gapped Hamiltonians given in Definition^ Then, there exists a linear operator 
, s , with the following properties: 

(i) F^^Hs) = H s and, hence, 7"y )0 (-ffo) = H 0) 

(ii) [Fji ,s(O u (r)), Po(s)} — 0, for any operator O u (r) supported on b u (r), u G A, 
(Hi) \\T Y AOu(r))\\ < \\O u (r)\\, 

(iv) Jy, s (O u (r)) = J2r'>o-Fy ,s(r';O u (r)), with J"y ;S (r'; O u (r)) supported on b u (r + r') and \\ Jy !S (r'; O u (r))\\ < 
\\O u (r)\\ fo (r'), for a rapidly decaying function /q . 

(v) 7" 7 /, s (Q„) - J" 7 ',o((3u) = J2 r >2 x ^\r), with Xu\r) supported on b u (r) and \\X^\r)\ < sJ\\Q u \\ fi(r), for a 
rapidly decaying function f^. 
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Proof. Define the operator J-y^ s as follows: 



T Y ^Ou(r))= / dtw Y (t)T t H '(O u (r)), 



(12) 



where t b (A) = e ltB Ae %tB and uy(t) is the filter function studied in [21 Q2], with the following properties: 

1. Compact Fourier Transform [22:: wy(u>) — 0, for |cj| > 7' and w) 7 '(0) = 1. 

2. Rapid decay: < w y {t) < 2e 2 i(y\t\) exp (~f ln 2 7 ( y| tD ) , for > e 1 /^ an d < w Y (t) < Y/tt, for all f <= 
Note that, we immediately get: 



FyAHs) = I dtw Y (t) e? tH >{H s )e- itH ° = ( f dtwy(t)\ H s = H„ 

J — OO \J — OO / 



(13) 



since J^^dtwy^t) — wy(0) — 1. Moreover, the new interactions commute with Pq(s), since for any ground state 
IV'o(s)) of H s and any eigenstate \i/j n (s)) orthogonal to Pq(s) with eigenvalue E n (s), we have: 

Jv,.(0«,(r)) |Vo(s)) = d>y ~ ^bOO) (^»( s )l °«( r ) IM*)> - 0, (14) 

where E n {s) — Eq(s) > 7', by assumption. This implies that: 

(I - P (a))T Y , a (O u {r))P (s) = = P (s)F r , s (O u (r))(l - P (s)), 

hence J r 7 / jS (0, i (r))Po(s) = Po(s)^ r 7 ', s (O u (r))Po( s ) = Po( s )-7 r 7 , ,s(O tl (r)), as desired. We also have by the triangle 
inequality: 



\\T Y , s (Ou(r))\\ < / dtw r (t) \\e^'(O u (r))e-^\\ = ||O u (r 



(15) 

Finally, the quasi-locality of Ty , s (O u (r)) follows from writing Jy )5 (0 u (r)) = J2 r , >0 Jy,^?"'; O u (r)), where for r' > 1: 



/OO 
dtuy (*) (rf \0 u {r)) - r t H ° {r+r - X) (O u (r))) , 



and 



7V, s (0;O„(r))= / *^(t)rf» (r) (O tt (r)) 



with the Hamiltonians in the exponent defined as follows, for < q < L: 



v:b v (l)Cb u (q) 



(16) 



(17) 



(18) 



Note that H™(L) = H s and that Jy, s (r'; O u (r)) has support on b u {r + r') and norm decaying rapidly in r' . The 
decay follows directly from the assumption on the almost-exponential decay of my(i) for \t\ large enough and, rapidly 
decaying Lieb- Robinson bounds in r' , on the norm of the commutator: 

[H-(r + r') -H»(r + r' l),rf ^ r+r '\o u (r))} , 
for small enough \t\. In particular, note that for 



we have: 



S s , t (r';O u {r)) 



S s , t (r';O u (r)) = 
1*1 



(Ou(r)) 



fl^(r+r'-l) 



(Ou(r)) 



_ff"(r+r'-l) / H"(r+r') 



(O tt (r))) - O u (r) 



< I \\[H"(r + r')-Hnr + r>-l),T^ r+r '\0 u (r))]\\dt> 



(19) 
(20) 
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Denning I q (u) = {Z C b u (q) : Z n (b u (q) \ b u (q — 1)) 7^ 0}, we continue the above estimate as follows: 



E / \\[Q v ,T«° [r+r '\O u {r))]\\dt' 



&„(i)ei r+r ,(u) 

+ E 

6„(fe)ei r+r ,(«) 



I'l 



[V v {k),T»° (r+r '\0 u {r))]\\dt>. 



(21) 
(22) 



Note that the norm of the commutator II \_Q V , t, 



H"(r+r') 



(0 M (r))] || can be bounded directly using Lieb- Robinson bounds 



for the evolution tJ 



, with rapid decay in the distance between the supports of b v {\) G I r+r i{u) and b u (r), which 



is at least r'— 4. Similarly, || [ V v {k),T^ s ^ r+r ' (0„(r))] || will decay rapidly in the distance between b v (k) e J r+r ./(u) and 

bu( r )j which is at least r' — 2k. Of course, for k > r'/2 we may use the simple bound || [V v (k), T ^ {r+r,) (O u (r))] || < 
2||K(fc)||||O u (r)|| < 2/ (fc) ||0„(r)||. On the other hand, we also have the trivial bound for all t E R: 



SsAr';O u (r)) < ||rf " (r+r,) (O u (r))|| + ||rf "^'"^ (O u (r)) || = 2||O u (r) 
Combining the above bounds with the usual technique of bounding the integral: 



(23) 



f(t)-g(t)dt<\\g\\ / \f(t)\dt+\\f\\ / \g{t)\dt 

-oo J-T J\t\>T 



\t\>T 

where |j/|| is the supremum of / over the appropriate domain, gives us the fastest decay for the optimal choice of the 
cut-off time T. In our case, f(t) = S s .t(r'; O u (r)) and g(t) = w 1 i(t). 

Now, for the final statement of the Lemma, we note that DuhamePs differentiation formula implies: 



/OO />S 
dt wy (t) / d s >r" B ' (Q u ) ds' = 
-co JO 

I 00 dtW Y (t) f if ( T^'{V)dt',T?°'{Q u ) 
J-co JO l J0 



ds' = 



J ds' J dt w Y (t) rf>' ( J i [t H °; (V) , Q u ] dt') , 



(24) 
(25) 
(26) 



where we performed a change of variables and used the linearity of the commutator to get the last line. Note that in 
the last line above, we have an evolution t^ s> applied to an integral involving the commutator: 

[r H p'(V),Q u } = ^[T H p'(V v (k)),Q u }. 

(v,k) 

It should be clear at this point that by applying the "telescoping sum" procedure, described earlier in this proof, to 
the evolution t H ^,' , each commutator in the above sum is itself a sum of terms supported on larger and larger sets, 
with decreasing norm. To see this, note that we may bound the norm of each such term trivially: 



|[' 



Hl,(k+k') 



(V v (k)) 



H\(k+k'-l) 



(V v (k)),Q u ]\\ <25 s ,,_ t ,{k';V v {k))\\Q v 



(27) 



using the appropriate bounds on 5 s >^ t ' {k'\ V v (k)), as described above. In all cases, we note that the commutator 
vanishes for k + k' < d(v,u) — 3, where d(v,u) is the distance between the sites u,v € A, so we only need to 
consider k' > max{0, d(u, v) — k — 2}. Finally, after localizing the terms inside the evolution t^ 3 ' , we perform another 
"telescoping sum" decomposition on this final evolution, arriving at the final decomposition of strength J (because of 
the norm of the terms V v (k)) and decay given by a function f[*\ whose decay depends implicitly on the decay of uy 
and Lieb- Robinson bounds corresponding to the unitary evolutions under Hamiltonians with interactions coming from 
Q u , u e A and V v (k), v e A, k > 0. The linear dependence of the bound on s follows from the triangle inequality: 

f ds'A(s') <s- sup POOH. 
Jo s'e[o,s] 



□ 
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B. Absorbing the perturbation. 

In this subsection, we focus on unitarily transforming the Hamiltonian H s — Hq + sV, satisfying the conditions 
of Definition |4j into a Hamiltonian H' s = H Q + V a + E s l, with the new perturbation Vg satisfying certain stronger 
conditions with respect to the local ground states of Hq, as discussed at the beginning of this section. Since the whole 
process is a unitary transformation combined with a global energy shift, the spectral gap of H s is identical to that of 
Hq + Vg. The transformation is done in four steps: 

1. We transform the terms Jy s := J" 7 ', S (Q U ) + sJ r 7 ', s (K) of H s = E«-7y, s > satisfying [Jy s , P (s)] = 0, into 
terms X u = U '(s) t ' T™, S U \s) - J" 7 ',o(Q«), satisfying: 

[X U ,P ] = U(s^[J^, rS ,P (s)]U(s) - [Jy,o(Ou), ft] = 0, 

where the unitary U(s) satisfies U (s) Poll (s) 11 = ft(s). At this point, the Hamiltonian H s has been unitarily 
transformed into H' s := U(s)^H a U(s) = P + Eu X u- 

2. We decompose X u = + xi 2) + x£ 3) , where: 

X^ = ^'AQu) ~ ^y,o(Qu), = U(s)^ 7 ,, s (Q u )U(s) - F 7 >,s(Qu), = sU{s)^ Y , s (V u )U{s), 

and show that each term is a strength J perturbation with decay given by functions /} , f^ and f[ , respec- 
tively, using Lemma [T] and properties of the unitary U(s) [3[ I17j. 

3. Using [X U ,P ] = 0, we rewrite X u = {l-P Q )X u (l-P )+A u +c{X u )l, where X u = X u -c{X u )l, A u = P X U P Q 
and c{X u ) = Tr(P X u )/Tr P . We show ||A U || decays rapidly in L* , using Local-TQO and the decomposition 
of X u from the previous step. 

4. For a (J, /i) perturbation X u satisfying [A„,P ] = 0, we prove that W u := (1 — Po)X u (l — P ) is a (J, w) 
perturbation, satisfying W u (r)P bu ^ = 0, with the decay w(r) closely related to fi(r). 

Note that throughout the transformation process, the perturbations X u , W u and A u depend implicitly on the param- 
eter s. Nevertheless, the locality and bounds on the strength of the perturbations are independent of s (some bounds 
have a linear dependence on s, but since s € [0,1], this makes the bounds only stronger.) 
We state the end result of the above process as a Proposition: 

Proposition 1. Let H s be the family of Hamiltonians considered in Definition^ Then, there exist unitaries U(s), 
such that: JJ(s) 1 ' H s U{s) = Hq + W + A + E (s)l, with W a (J, w) -perturbation satisfying, W u (r)Pb u M = 0, ||A|| 
rapidly decaying in L and Eq(s) = Tr(P s Po(s))/ Tr Pq(s) . 

Proof. Using Lemma [T] we see that the first step implies H' s := U(s)^ H s ll(s) — U(s)^J r 1 ',s(H s )U(s) = J- y ,q(Hq) + 
E„ X u = Hq + E u X~u, where [X u , Po] = 0. The second and third step imply: 

H' S = H + 2(1 - Po)X u (l — P ) + A + c(A)l, 

u 

for A = E«A U a vanishingly small perturbation and X = ^2 U X U . Noting that X = U(sy H s U(s) — Hq and 
Po{s) = U(s)P U(sY , the cyclicity of the trace implies: 

c(X) = Tr(if s f/( S )P f/( S )t)/TrP - Tr(P P )/ Tr P = Tr (P s P ( S ))/TrP ( S ), 

so H' s = Hq + (1 — P )A(1 — P ) + A + E (s)t. Finally, setting W = E u W u and using the second and fourth steps, 
we get the desired result. □ 

We focus, now, on proving the last three steps above, since the first step follows immediately from the properties of 
.F 7 ' )S and the quasi-adiabatic evolution U(s) (see, e.g. [21 Prop. 2.4].) For the second step above, we make extensive 
use of results obtained in [3[ [17] . We highly recommend reading through [T7j to get a clear conceptual understanding 
of the ideas behind Hastings' quasi-adiabatic evolution U(s) and its many applications to problems in condensed 
matter physics and quantum information theory. Below, we reference results from [3], from which one may extract 

explicit bounds on the decay of X^ and X^\ The properties of Xu^ follow directly from Lemma [l] 



10 



Lemma 2. Let H s be the family of gapped Hamiltonians described in Definition^ and define Gy, s (O u {r)) '■= 
a^(O u {r)), where is defined in £3j/ as the spectral flow corresponding to H s and O u (r) is an operator supported on 
b u (r). Then, 

L-r 

g Y . s (O u (r))= ]T £V, s (r';O u (r)), r > 0, (28) 

r'=r 

with Qj> lS (r'; O u {r)) supported on b u (r + r') and 

\\g Y Ar';O u (r))\\<J\\O u (r)\\g(r r ), r' > r , (29) 
||0 7 ', s (ro;O»(r))-O„(r)|| < sJ \\O u (r)\\ g(r ), ||0y )S (r o ; O u (r))\\ = \\O u (r)\\, 

for a rapidly decaying function g. 

Proof. We introduce the decomposition Gy , s (O u {r)) — X)r'=r Gy, s (r'; O u {r)): 

G Y , s (r / ;O u (r))=a^(O u (r))-as r '- 1 (O u (r)), r' > r + 1, (30) 
GyAro;O u {r)) = as r °(O u (r)) r > 0, (31) 

where we have defined A„ = b u (r + n), n > and used Eq. (4.42-4.46) of [3] to define aJ(A) = U Y (s) A U Y (s): 

d s U Y (s)=iD Y (s)U Y (s), U Y (0) = 1, D Y (s) = J2 *a(Z,s). (32) 

ZcY 

First, notice that since each term ^\(Z, s) is supported on Z by definition (see [3]), the terms G^', s (r';O u (r)) have 

support on b u (r + r'), as desired. Moreover, since Gf',s(fo;O u (r)) — a s r °(O u (r)) = (s) O u (r)U^ ro {s), we 

immediately get ||£/ 7 ', s (ro; 0„(r))|| = ||0„(r)||. To bound the norm of Qy s (r'; O u {r)) for r' > r , we need to bound 

\\ a s r ' {O u (r)) — a^ r ' _1 (O u (r)) || . As we show in Appendix [aJ the norm \\a^' J (O u (r)) — a^ r ' ~ 1 (O u (r))\\ is bounded 
above by: 

J\\O u (r)\\g(r>), (33) 
for a rapidly decaying function g. A similar estimate shows that: 

\\Gy, s (r ;O u (r))-O v (r)\\< f ds' \\[D Ar J S '),at (O u (r))}\\ < sJ \\O u (r)\\ g(r ). (34) 

Jo 

□ 

The third step makes use of the Local- TQO condition, allowing us to extract the groundstate energy of H s as a 
single quantity Eq(s), up to a small correction. In other words, the following Lemma asserts the vanishingly small 
splitting of the groundstate subspace for the gapped Hamiltonians H s . 

Lemma 3. Let X u be a (J, fi) -perturbation centered on u G A. Then, Assumption [7] implies for A u = PqX u Pq, 
X U = X U - c{X u )l, c(X u ) = Tr(X u P )/TrP : 

L" L 

\\A u \\<j(j2fi(r)&o(L-r) + 2 £ Mr)). (35) 

r=0 r>L' 

Proof. By assumption, X u = J2r>o X u (r), with X u (r) supported on b u (r) and HA^r)! < J f\{r). We introduce an 
energy shift by redefining the local terms: 

X u {r) := X u {r) - c(X u (r))t, c{X u (r)) = Tv(P X u (r))/ Tr P . (36) 

From Local-TQO we know that for r < L* , 

\\P X u (r)P \\ < \\X u (r)\\ A (L - r). (37) 
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On the other hand, for r > L* , we have: 

\\P X u {r)P \\ < ||^«(r)|| + \c(X u {r))\ < 2J/i(r). (38) 



Using (37 1 and (38), we get the desired bound: 

L" L 

\\PoX u P \\< j(X)/i(r)Ao(£-r)+2 ]T A(r)). (39) 

r=0 r>L* 

□ 

We now move to the last step of our transformation. 

Lemma 4. Consider a frustration-free Hamiltonian Hq and a perturbation X u with strength J and decay fi(r), 
satisfying [X u ,Pq] = 0. Set W u = (1 — Pq)X u (1 — Pq), where X u = X u — c(X u )l. Then, W u is a strength J 
perturbation with decay w(r) given below and W u (r)P^ r ^ = 0. 

Proof. By definition, W u Po = 0. Moreover, using [X u , P ] = 0, we have: X u = W u + A u , where A u = P X U P . This 
implies that: 

L-l 



W u = J2 X u(r) + (X u {L)-A u ). (40) 

r>0 

Let us choose I < L* . We apply Corollary [2] to J2t=o X ^ ( r ) : 

e it 
\\J2Xu(r)P K ( 2 e)\\ < ||^l„(r)Po|| +2||^l„(r)|| y/E^jtj (41) 

r— r— r— 

L I 

< \\W U P \\ + || Y,X u (r)P \\ + || A„|| + 2 || Xu(r)\\ y/A (e) (42) 

L £ 

<2j(^/ 1 (r) + 2^/ 1 (r) v / A^))+||A„|| ) (43) 

where we used ||X u (r)|| < 2||X u (r)||. We turn our attention to the final decomposition of W u into interaction terms 
that are annihilated by local ground states. 

Define an orthogonal unity decomposition 1 = X)m=\ E m , as m : 

Em = Pb u (m-1) ~ Pb u (m), 1 < m < L + 1, P u = 1, Pb u (L+l) = 0. 

Using T,i< P ,r<L = Ej 2 > 2 E P+ r=j and E r <„< s E p = A„(r-i) - A„( s ), we get: 

L-l L 

W u := (1 - P )X U (1 -P ) = J2 YuU) + Y U (L) + ]T Z. u (2g - 1) + Z U (L), 

3 >2 2q=2 



where we defined: 



Y u (j) := E v E ^> 2 < j < L — 1, 

p+r— j q—0 

Y U (L):= J2 E P Y J X u {q)E r , g(p,r) := Lmax(p, r)/2j , 

p+r— L q—0 

Z u (2g - 1) := (1 - P K{2q -i)) X u (q) (1 - P M2g -i)), 2 < 2<j < L, 
Z tt (L) := (1 - P ) £ X u (q)(l-P ). 

q>L/2 
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Note that Y u (j)P b .( u ) = and Z u (2q — l)Pb u (2q-i) = 0. In other words, Y u (j) and Z u (2q — 1) are supported on 
bj(u) and b u (2q — 1), respectively, and are annihilated by the local ground states. The same holds for the terms Y U (L) 
and Z U (L). Finally: 



\\Z u {2q- 1)|| < \\X u {q)\\ < 2Jf 1 (q), \\Z U (L)\\ < 2 J £ / x (g), 

q>L/2 

£(p,r) 

\\Y u (j)\\ < \\ E p E Xu{q)E r 

p-\-r=j q—0 

LP/2J Lp/2J 
^ 2 E E ^ 2 E \\ P K( P -l) E 

P>j-P <?=0 P>j—P 9=0 

3-1 i LP/2J 

^ 2J E ( E /iW+ 2 E /iWVMrfj+iiiA, 

p>Li/2j r>[_p/2\ r=0 



(44) 
(45) 
(46) 
(47) 



where the last line follows from (43). To bound the norm of Y U (L), we simply sum over the bounds for the norms of 
Yu(j): 



2L j-1 



LP/2J 



|y„(L)||<2j]T Yl ( E AM + 2 ^ A(r)yA^))+2L 2 ||A u | 

3>Lp>[j/2\ r>[_p/2\ 



r=0 



We define the function w(r) as follows: 



w(r) = \\Y u (r)\\, r = (mod 2), r < L, 

«j(r) = ||y u (r)|| + ||2' u (r)||, r=l (mod 2) r < L, 

w(L) = \\Y U (L)\\ + \\Z U (L)\\. 



(48) 



(49) 



For Ao(r) and fi(r) decaying rapidly, we see from ( 44p8| and Lemma [3] that w(r) is also rapidly decaying. □ 



V. PROOF OF RELATIVE BOUNDEDNESS 



We now generalize the proof of Proposition 1 in [6] . 

Proposition 2. Let W be a (J, w) -perturbation satisfying: T / F tl (r)Pf, u (r) — 0. For Hq a frustration-free, locally-gapped 
Hamiltonian on A C 1> d , with gap decay 7(7*), there exists a constant 

c:=C, dE4H E ^0)' = r- o <n < ••• <r M = L, (50) 
^ — ' 7(rfc) ^ — ' 

where Cd is the volume of the unit ball in Z d , such that for arbitrary states \ip): 

\ (ip\W \ip) \ < c • J (ip\ H Q \ip) . (51) 

Proof. First, we partition the lattice A into disjoint sets Cj(r), I < j < N(r), such that any two distinct points 
u,v E Cj(r) satisfy b u (r) n b v (r) = 0. We treat the case 8r > L separately, using the trivial decomposition Cj(r) = 
{uj}, Uj <E A, 1 < j < N(r) = L d . To accomplish the decomposition for 8r < L, we start by dividing A into an even 
number of boxes of size 2r', in each direction, such that 2r < 2r' < 3r. In particular, we have L — 2rq+ro, < ro < 2r, 
so we may choose 2r' = 2rq/(q — b) + ro/(q — b), where b = q (mod 2). Since we assumed 8r < L, we know that 
2r • 3 + ro < L => q > 4, which implies the bounds on 2r' . Note that 

L/2r =q-b = (mod 2) 

and that r need not be an integer. 

We build the partition one dimension at a time, as follows: In one dimension, we "color" every other interval of 
length 2r' on the line with the same color. So, in one dimension, we have 2 colors. In 2 dimensions, we fill the plane 
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FIG. 1. Partition of the lattice A into an even number of boxes of size 2r' each. The black dots indicate the points Uj inside 
some set Cj(r). The label a corresponds to the label of the box in which the point is found. 



with an even number of rows of 2r' x 2r' squares, where every other row looks the same (same "color" pattern) and 
within each row we have an even number of alternating colored boxes (see Fig. [T]) We can continue like this for any 
dimension d of our lattice, requiring 2 d colors in total. Now, because of the way we have constructed our partition, 
any two points on distinct boxes of the same color will be a distance at least 2r' > 2r apart, as desired. By choosing 
exactly one point Uj from each box B a of the same color, we get a family of disjoint sets Cj(r) that cover A, each 
satisfying the desired property of internal sparsity. The total number N(r) of these sets is bounded by the number of 
colors multiplied by the volume of each box. So, we have: 



N(r) < 2 d {2r') d < 6 d r d 



and 



N(r) 



A = (J Cj(r), with Cj(r) = {u", a > 1} defined above. 

4>i 



From this point onward, we set: 



bj : a( r ) : = b u <>{r), Uj G Cj(r), a 



> 1. 



(52) 



(53) 



(54) 



Note that [P^ a ( r ), P^. 6 ( r )] = 0, since the projectors have disjoint supports for a ^= b. Define an orthogonal decompo- 
sition of unity as follows: 



IT (y^-Pb jia (r)) + (l-ya)P bi ^r)), 3^ S {0, 1} 



It is easy to check that Y,y ' W = 1 and Rf( r )Rf( r ) = 5y,zRj(r). 



Before continuing, wc define a decomposition of the perturbation W as follows 

W 



M N(r k ) 

E E ^( r ')> Wj(r k ):= Yl W u a(r k ), 

k>l j>\ a:«°6C,(r t ) 



W u {r k ):= W u{r), = r < n < r 2 < . . . < r M = L. 



(55) 
(56) 



r>r k -i 
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Since W(r) is locally block-diagonal, we have 

W u? (r)P bj:a( r) = W u «(r k )P bjArk) =0. 
Moreover, recalling that w(r) := sup u ||W tt (r)||/J, we have: 

11^.(^)11 < j Yl ( 57 ) 

r>Tfe_i 

Setting w(r k ) ■= X)r>r fc _i w ( r )i we wm show that: 

| (V| T^-(r fe ) |V) | < J - u>(r fc ) (V| G^r*) |V) , (58) 

where, 

G,(r fe ) :=]>>• E ( 59 ) 
n>l \y\=n 

Define Wj^ a {r k ) ■= W u <*(rh), for it| G Cj(r k ). By definition of Ry(r) and the local block-diagonality of ^(rfc) = 

Ea^>( r fc), we have: 

^(r fe )# 7 (r fc ) J Rf(r )fc ) = 5 3 ;,2 ^ J^(r fc )^, a (r fc )i^(r fc ) => (60) 

\\Rf(r k )W 3 (r k )Rf (r k )\\ < Jw(r k ) 6y tZ ■ \y\, (61) 
where the equality in the first line follows from the observations: 

1. W j , a (r k )Rj(r k )=0, for y a = 0. 

2. If y ^ Z, then either Z a = when y a = 1 and W^r^-Rf (r fc ) = 0, or y b = and Z h = 1, for b ^ a, so 
R? {r k )Wj, a (r k )Rf (r k ) equals: 



The bound (61) implies for arbitrary states 

i (vi w^(r fc ) iv) i < E i m ^w^w^^) w I 

< E ll^(^)^(r fc )^(r)|| (VI i#(r fc ) |V) < -M(r fc ) £ M M ^(r*) IV) 

y y 
= Jw(rk){i>\G j (r k )\ip). 

Recalling that H b . a{r) = E« : 6„(i) C 6 J , (l (r) denne H j( r ) = Y* a H b iia {r) and note tnat: 

JV(r) 

^ H s (t) = (E M W 1 ) c u a>i6j,a(r)))Q« < Q/ffo, (62) 

3=1 «6A 3>1 

where is the volume of the unit ball in 7L d and Ind denotes the indicator function. To see this, note that 
5Zj>i Ind(6 u (l) C U a >ifej. a (r)) is bounded above by the number of balls of radius r that contain b u (l), since 
Ind(fe„(l) C U a >i&j. a (r)) is non-zero only when j is such that bj^ a {r) := b u ^(r) includes b u (i) for some unique a, 
since bj ya {r) n bj_ b (r) = 0, for a =/= b. Thus, the indicator sum is bounded by the number of places we can put 
b u (l) inside a ball of radius r, which is the same as the number of balls of radius r containing b u (l). The former is 
upper-bounded by the number of places we can put bo(u) = u inside a ball of radius r, which is exactly equal to Cd r . 
Using the assumption that Hq is locally gapped and frustration-free, we have 

H bja[r) > 73>(r)(l - P bj , a (r)), y a = H b] a[r) R y {r) = 0. 
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Then, we can proceed as follows, for arbitrary states \ip): 

&\H j (r)\i>)=J2 E W R *(r)H hjAr) R?{r)H) 
y a-.y a =i 

>E E iar)mRj(r)(l-P bj . a (r))Rf(r)W 
y a-.y a =i 

< 7i(r) E n E M R K r ) W =► Gj{rk) < H^hiiru), 

n>l \y\=n 

where 7j(r) := min a 7j, a (r). We have just shown that: 

| y>\ Wjin) |V) | < Jw(r k ) (M Gj(r k ) |V) < JwinJ/jjin) (VI fl^r*) |V> . 

This implies that: 

M N(r k ) M N(r k ) 

kviwiv)i<E E iw^wi^i < jE*w E (vi^(^)iv)/7,(r fc )- (63) 

k=l j = l k=l j = l 

Setting j(r) := min.,- 7j(r), we have: 

\(il>\W\il>)\<cJ(fl>\H \fl>), (64) 



where c = C d J2k=i r k^( r k)/li r k), using (52) and (62). □ 



We note a few things about the definition of the structural constant c appearing in Proposition [2] (i) The choice 
of the rk : < r% < . . . < tm — L is left unspecified beyond the boundary conditions, so that one may optimize the 
constant c given the relevant decays w(r) and 7(r). For example, if j(r) decays particularly fast for some values of r, 
we may want to avoid these values by summing over them within the rapid decay of w(r). (ii) Assuming polynomial 
decay of 7(7^) for the selected values of r*., it is clear that for fast enough decay of w(r) (superpolynomial, or even 
fast enough polynomial decay), c is a constant depending only on the dimension d of the lattice and the decay rate of 
7(r&) and w(r). 

VI. FINAL ARGUMENT: 

The argument below follows the bootstrapping argument in [B], putting together all the pieces to prove stability 
of frustration-free systems under weak perturbations. 

Theorem 1. Let Hq be a frustration-free Hamiltonian, with spectral gap 7, defined in Section\H\ satisfying the Local- 
TQO and Local-Gap conditions. For a (J, fo) -perturbation V, with /q rapidly- decaying, there exists a constant c, 
given in PropositioniOA such that for J < ^- and L large enough, the spectral gap of H + V is bounded from below by 
7/2. 

Proof. Let Jo = 7^0 where c is defined in Proposition [2] Assume that for J < J , H + sV has gap 7 S > 7/2, for 
< s < So < 1 and 7 S < 7/2, for sq < s < s\ < 1. Then, we have proven in the previous sections that for < s < sq: 

U{s)\H a + sV)U(s) =H + W s + A s + E s -t, 

where U(s) generates the spectral flow U(s)P (0)U(s)' < — Pq(s), and | (ip\ W s \ip) \ < c- J (ip\ H \ip). Moreover, E s ■ 1 
is an overall energy shift and the norm of A s decays rapidly in L. We note at this point that the Local-TQO cut-off 
L* plays a crucial role in the decay estimates of ||A S ||, for example. In particular, L* should scale as L a , for some 
a G (0, 1], in order for estimates like the one in Lemma|3]to decay rapidly in L, especially since the estimate for ||A S || 
is given by L d su PueA ||A„|| - L d max{A (L*), fo{L*)f. 

We will show that, for large enough L, the gap of H + W So + A So , which is 7 Sf) , is strictly greater than 7/2, 
contradicting the assumption that 7^ < 7/2, for so < s < Si, implying that sq = 1. 

Now, we note that Ho + W s + A s = U(sy(Ho + sV — E s ■ l)U(s) has eigenstates in the subspace Po(0), with energy 
within ||A S || from 0. Those states correspond to the "groundstate" subspace -Po( s ) of the unrotated version of the 
Hamiltonian, Hq + sV. In general, eigenstates of Hq + sV in the Pq(s) subspace will have different energies, but only 



16 



up to vanishingly small error in the system size, L. To see this, note that for |^o(s)) = Po(s) \^q(s)) an eigenstate of 
H + sV with energy E (s), setting |#„) := U(s)^ |* («)), we have P (0) |*o) = U(s^P (s) |*„(s)) = |* ) and 

{H + W S + A s ) |* ) = (*)(#(> + *V - £ s • l)U(s) |* > = - E.) |*„) • 

Moreover, since | (* | W s |* ) I < cJ (* | #o |*o) = 0, the energy £b(s), satisfies 

(flo + W s + A s ) |¥„) = A s |* > = E Q (s) |* ) . 

Hence, |i?Q (s) | < ||A S || < 7/13, for large enough L. 

Moreover, any state \ip±) orthogonal to the Pq(0) subspace has energy at least (1 — cJ)j— ||A S ||, which follows from 
the linear relative bound on W s and the gap of H : 

(<fa\H a + W s + A s |^) > (falHo - | (ih\W, | - | (ipi\ A * I 

> (1 - cJ) H - ||A S || > (1 - cJ) 7 - ||A S ||. 

Now, taking the linear size L large enough (such that ||A S || < 7/13) and recalling that cJ < 1/3, we see that the gap 
is bounded below by: 

7s > (1 - cJ) 7 - 2||A S || > (20/39)7 > 7/2, 

for all s S [0, so]. But, by definition of so, we know that for any < e < s% — sq, we have j(sq + e) < 7/2, a 
contradiction of the continuity of 7 S for finite L. Hence, So = 1 and 7(1) > 7/2, for all J < Jo and L large enough to 
satisfy ||A S || < 7/13. 

Finally, this implies that Hq + sV has gap 7 S > 7/2, for all s E [0, 1] for the parameters of Jo and L given above 
and determined implicitly by the decay fo(r) of V, the Local-TQO decay Ao(r), the local gap 7(7-) and the linear size 
L of the system. □ 

VII. DISCUSSION 

We have shown that topological phases corresponding to the groundstate sector of gapped, frustration- free Hamilto- 
nians are stable under quasi-local perturbations. In particular, the splitting in the energy of the groundstate subspace 
is bounded by a rapidly-decaying function of the system size and the gap to the high-energy sectors remains open for 
perturbation strengths independent of the system size. Our result also implies the stability of symmetry-protected 
sectors under perturbations that obey the symmetry of that sector. To see this, note that the Local-TQO condition 
is only required for operators Oa obeying the symmetries under consideration. Stability follows using the methods 
already found here, as is already discussed in [6]. 

Moreover, it can be shown [3S] that parent Hamiltonians of Matrix Product States (MPS) satisfy the Local-TQO 
condition with Ao decaying exponentially in these cases. Combined with the results of Nachtergaele |28j and Spitzer et 
al. |42j on the spectral gap of one-dimensional, frustration- free Hamiltonians, our result shows that parent Hamiltonians 
of MPS have stable low-energy spectrum against arbitrary, sufficiently weak, quasi-local perturbations. 

At this point, there are several important open questions that one may want to pursue: 

1. First, can we extend the one-dimensional methods, introduced by Nachtergaele in [28], to prove rigorous spectral 
gaps for frustration-free systems in two dimensions? 

2. Can we use the decay rate Ao in Local-TQO to classify different phases of matter? In other words, is there 
a canonical parent Hamiltonian construction, given a groundstate subspace, such that Local-TQO is satisfied 
with an optimal decay Ao? For example, see [51 Sec. 2.4], where two different Hamiltonians for the toric code 
give very different decays Ao; one being constant and the other being a step function. 

3. Under what conditions (maybe translation invariance, and/or scale invariance) is Local-TQO satisfied, if the 
global groundstate space is topologically ordered? This is related to the previous question. 

4. Can one rigorously prove that all gapped, frustration-free systems with non-vanishing, local interactions, satisfy 
the Local-Gap condition, using the detectability lemma, or some other method? 

5. What features of parent Hamiltonians of PEPS can one use to prove that the Local-TQO and Local- Gap con- 
ditions are satisfied? 

6. Can one use the present result to prove the Haldane conjecture, by identifying a gapped, frustration-free Hamil- 
tonian which satisfies the Local-TQO and Local-Gap conditions, sufficiently close to the Haldane phase [231 ? 
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7. Can we identify two-dimensional, frustration-free systems satisfying Local-TQO that exhibit Topological Order 
at non-zero temperatures? What about three-dimensional frustration-free systems? 

8. Finally, what is the form of Local-TQO that is truly necessary (and sufficient) for stability of the spectral gap 
against local perturbations, for frustrated, gapped Hamiltonians? Is the Local-TQO condition always true for 
generic gapped Hamiltonians? 
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Appendix A: Lieb- Robinson Estimates for Lemma [2] 

We follow the proof of Thm. 5.2 in [3] to bound \\a b s uir ' ] (O u {r)) - a s u{r ' ~ 1] {O u {r))\\ in Lemma |J To use the 
notation of [3], we set A„ = b u (n + r), for n > 0, so that O u (r) is supported on the set Aq. We have: 



\\afr'{0»(r))-a£ r '- 1 (0 1l {r))\\ < / \\[D K ,(t) - D^tla^'- 1 (O u (r))]\\dt (Al) 

Jo 

<(2C f \\*\\f f g(t)dt)\\O u (r)\\ ]T £ F»(d(x,z)) (A2) 

+ 2 (2C f \\*\\f f 9it)dt\ \\O u (r)\\ ]T J2 F *(d(x,z)) (A3) 

+ 4«||*|| J ,||O u (r)|| Yl E Md(x,z)) (A4) 

zeA r'-l\ A L2(i-'-l)/3J xeA L('-'-l)/3J 

+ 16s ||0*||„F„(O) Y, VGi(n~l) \\O u (r)\\ |A L2(f .,_ 1)/SJ WK{\[r' - 1)/3J), (A5) 



n>l 

where the following functions are defined: 



with F(r) = (1 + r y {d+1) and, 



Mr)=u»(£-r)F(£-r), (A6) 



u^{x) otherwise. 



x 



In 



x 



u^(x) = cxp < — fj, — 2 — f for /i > 0, x > 1. (A8) 



Moreover, 



and 



2v a J av a 



K(x) < 4G 2 



(m + ^Pe^, (A10) 
V 2v a J a z vi 
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where C?2(C) is defined for ( > by 

1 f7,354 <C< 36,058 

2(0 ~ f ' \l30e 2 C 10 u 2/7 (C) C> 36,058. " 

Finally, the remaining constants (e.g. Cp = 2 d+1 J2rei d ^( r )' ll^ll-F' II^IUj ^a(0)i a, v a , etc.) are defined in Sec. 
4], where one also sees that: 

2011*111?/ g(t)dt < e 25 "*"^ - l. (All) 
Jo 

We note that the decays involved (i.e. the functions K(x) and F$(r)) are most likely not optimal, but we consider 
it an important exercise to give a rough estimate of the order of magnitude and the dependence to the spectral gap 
7', appearing in Definition |4j of the constants involved in the decay of the perturbations we are studying. In any 
case, the relevant decay for our computation in Lemma ^ is dominated by tail estimates of the function u^{r'), which 
defines a sub-exponential decay with [i ~ 7'. Here, we nave sketched how one may derive rigorous estimates, using 
the tightest Lieb-Robinson bounds available at this point in time. 
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